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Abstract 

In this paper we present the first steps for obtaining a discrete Quan- 
tum Mechanics making use of the Umbral Calculus. The idea is to dis- 
cretize the continuous Schrodinger equation substituting the continuous 
derivatives by discrete ones and the space-time continuous variables by 
well determined operators that verify some Umbral Calculus conditions. 
In this way we assure that some properties of integrability and symme- 
tries of the continuous equation are preserved and also the solutions of the 
continuous case can be recovered discretized in a simple way. The case 
of the Schrodinger equation with a potential depending only in the space 
variable is discussed. 

1 Introduction 

In the last years the interest in physics for discrete space-time theories has 
increased mainly in relation with quantum gravity ^IJ. In a scenario where the 
physical space-time may be discrete (i.e. there is a fundamental length in the 
space-time -let a be the fundamental space length and r the fundamental time 
length - that could be related with Planck's constant), continuous theories 
would only be approximations to the reality. Hence, it seems interesting the 
study of discrete physical theories such that in the limit of a going to zero we 
recover the well known and established continuous physical theories. One way 
to do that is to discretize continuous physical theories, in particular quantum 
theories. It is obvious, that after discretization some properties will be conserved 
and other ones will disappear. A detailed study about when, how and why this 
happens would be pertinent and very interesting. This is the aim of this work: 
to study how to discretize a physical theory in such a way that properties related 
with symmetries in the continuous case can be preserved in the process, and so, 
to analyse the behaviour of quantum mechanics after discretization. 

The Umbral Calculus is the mathematical tool that we will use to discretize 
Quantum Mechanics. It has been used recently to provide discrete representa- 
tions of canonical commutation relations (i.e. [9a;,a;] — 1) in order to discretize 
linear differential equations in such a way that the continuous point 

symmetries are preserved the most cases. 

In this paper we will consider a fixed lattice and use the umbral correspon- 
dence to obtain discrete solutions maintaining certain commutation relations 
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The history of the Umbral Calculus goes back to the 17th century: Taylor, 
Newton, Barrow, Vandermonde, etc. jlOj found some theorems and expressions 
relating powers with polynomials. However, in the second half of the 19th 
century through Sylvester, Cayley and Bhssard appeared the terms 

'umbrae' (shadow in Latin) and 'Umbral Calculus' for the first time in relation 
with a set of 'magic' rules of lowering and raising indices. Some mathematical 
developments, like the theory of Sheffer polynomials (called Appell polynomials 
or more generally polynomials of type zero) |14|-[17| and the theory of abstract 
linear operators by Pincherle [13], converged to Umbral Calculus. Finally, in 
the second half of the 20th century Rota, Roman and collaborators jl8j-[23] 
developed the Umbral Calculus as a linear algebra of operators. Thus, they 
formulated the theory of the umbral algebra and its dual in the space of formal 
series. A large bibliography about Umbral Calculus and its applications can be 
found in Ref. [Mj. 

About recent results, we can mention that new methods in Umbral Calcu- 
lus have been introduced for computation of series by integration [25 . Also 
Ref. |26|-[30j study the symmetries of nonlinear difference equations or multi- 
dimensional ones related with Quantum Mechanics. 

This paper is organised as follows. In the next section we introduce the 
basic concepts of Umbral Calculus in terms of the theory of linear operators. 
In section |3] we show the main facts relative to the discretization of the (con- 
tinuous) Quantum Mechanics. The most usual representations for the discrete 
derivative operator (right, left and symmetric discrete derivatives) are analysed 
in section U] In order to see how the above umbral correspondences work we 
present in section [S] their action on some non-polynomial functions: the expo- 
nential and the trigonometric functions. The next section is devoted to present 
the basis for constructing a discrete Quantum Mechanics founded on the Um- 
bral Calculus. Thus we present an umbral version of the Schrodinger equation 
for some potentials. Finally some conclusions close the paper. 

2 Introduction to Umbral Calculus 

This introductory exposition of the modern theory of the Umbral Calculus, 
developed mainly by Rota and Roman, follows mainly Ref. [23]. The umbral 
algebra, £, is the algebra of linear operators acting on the algebra, J^, of formal 
power series in a variable or acting on the algebra of real or complex polynomials 
in a variable, V. The elements of J- are of the form 

oo 

/(x) = ^a„a;", a„ e F, 

ri=0 

and those of V are 

n 

Pnix) = ^ aix\ e F, n e N, 

being F a field of characteristic zero (in our case it will be identified with the 
real numbers M, but it could be, for example, the complex numbers). 
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The linear operators of £ called shift operator, T, and coordinate operator, 
X, act on the elements of J- or V as follows 

Tcr ^p{x + a), 

(1) 

X ■ p{x) — X p{x). 

All the results of this work can be extended to separable equations of any number 
of variables. For non separable equations the method is the same but the results 
may be quite different. 

An operator O S £ is said shift invariant if and only if commutes with the 
shift operator, i.e., 

[0,T,]=OT,-T,O^Q, VaeF. 
The so-called delta operators, A, are those operators of L such as 
A-a; = C7^0, ceF, 

(2) 

[A, T,] = AT, - T^A = 0, VfTGF, 

where c G F is any fixed constant. It is easy to prove that 

A • a = 0, V a e F. 

For every delta operator A there is a series of basic polynomials 
{Pq{x),pi{x),P2{x), . . . ,Pn{x)} of V verifying 

Po{x) = 1, 

p„(0)-0, Vn>0, (3) 

A -Pnix) = n Pn^i{x). 

We see that A acts as a lowering (index) operator on its associated polynomial. 
A very important fact of the theory is that each pair {A,pn{x)} defines a unique 
realization or representation of the umbral algebra. There is a bijective corre- 
spondence 3, between basic series of polynomials and their A operators. Hence, 
-3 may be defined with only one element of the pair. 

The commutator of an operator O with the coordinate operator X 

O' ^[0,X]^OX - XO (4) 

is called the Pincherle derivative of O. 

Given a delta operator A we can find an associated operator ^ such that 

[A,e] = l. (5) 

This operator ^ is not shift invariant and together with A and 1 close the 
Heisenberg-Weyl algebra or, in other words, expression ([S]) is a Heisenberg-like 
relation between A and its conjugate It is worthy to note that ^ is associated 
with the corresponding basic series of the operator A by 

r-i=Pn(x). 
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Although only in a determined correspondence • 1 coincides with the Pock- 
hammer symbol a;^"-*, from now on we will use the following notation 

C ■l=Pn{x) = X^^\ (6) 

which is easy to prove taking into account that 

Obviously, there is only one way to define ^ if its associated polynomial series 
obeys all the relations ([3|). Given an operator A its associated operator ^ must 
be 

where /3~^ is the first Pincherle derivative of A. It is easy to demonstrate that 
is invertible and (3 is shift invariant. 
Let Ai and A2 be two delta operators with associated ^-operators ^1 and 
^2, respectively. A map 7?. is called an umbral correspondence if 

7^ : ?r - C2 , Vn e N. 

The umbral correspondence also induces a correspondence TZ between two delta 
operators. Summarising, the following diagram is commutative 





n 


?2 






n 


Ai . 


n 


A2 












.1 = 4"^ 



From the above diagram one easily sees that there is a triplet {A, <^,p„(a;)} 
that defines a unique representation of the umbral algebra. Moreover through 
3 we can go in a bijective way from one element to another of the triplet in 
such a way that only one of these elements is necessary and sufficient to define 
the umbral representation. The bijective application TZ bring us from a umbral 
representation to another one. 

3 Umbral correspondence and discretized quan- 
tum mechanics 

As we mentioned before, the main goal of this work is the application of the um- 
bral correspondence to discretize Quantum Mechanics. In fact, the 1-dimensional 
Quantum Mechanics is an umbral realization via the identification 

A = 9„ ^ = X, p„(x)=x". 

Since we are interested in the discretization of Quantum Mechanics we will use 
the umbral correspondence to relate the triplet {dx, X, x") with another one 
defined on the lattice field = {x G M |a; = ma, m £ Z} with lattice parameter 
cr e M. 
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Let us start defining a general difference operator^ A, as a delta operator 
such that in the Hmit when a goes to zero becomes the differential operator dx ■ 
Supposing that A is a linear difference operator then it will have a polynomial 
dependence on the shift operator T, i.e., 

A = ^ E anT", (7) 

where iV G N, a„ G M and cr G R is the lattice parameter. Imposing condition 
(2^), i.e. A ■ X = 1, since condition ^p) is trivially verified we get 

EL-j a„ = 0, 
E„=-j nan = N. 

These conditions ^ guarantee that in the limit when <t goes to zero A goes to 
dx as it is easy to prove. 

Recalling that the pair of Quantum Mechanics operators {dx, X) verifies the 
Heisenberg relation 

[dx,X]^l, 

we call discrete position operator to a non-shift invariant operator ^ such that 
with the difference operator A verifies a similar Heisenberg relation, i.e., 

[A,e]-1. 

The umbral correspondence allows us to relate continuous linear equations 
h and their solutions / if they can be developed in power series, with umbral 
discrete versions h and /, respectively, as follows 

h{d^f{x),xn^O ^ /i(A"/(0,r) -1 = 0, n,sGN; 

CSO OO (9) 

This is the main result of Umbral Calculus that we will exploit here for ours 
purposes. Umbral Calculus allows us not only to discretize continuous differ- 
ential equations obtaining discrete equations but also to get solutions of these 
equations if their continuous counterparts are polynomial solutions or may be 
developed in power series which converge after the correspondence. Note, how- 
ever, that with the umbral correspondence we obtain an operator which must 
be applied to 1 to get a difference equation or its solution. 

For the particular case of a symmetric difference operator, i.e., a difference 
operator A„ such that 

A<,./(x) = -A_,./(a;), 
we have that the general expression ([7| reduces to 

^- = ]^E«"(^"-^"")- (10) 

n—l 
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4 Differential operators for discretized quantum 
mechanics 



In the following we will use some particular expressions (or representations) for 
the differential operator A, all of them of low (first and second) order on the 
shift operator T. Let us consider the two possible representations for A of first 
order, the right discrete derivative A4. and the left discrete derivative A_, and 
the (only) symmetric discrete derivative of second order A^. 

The right- correspondence is determined by the couple of operators (A+,^_|_) 
given by 

A+:=i(T-l), f3+^T-\ 
a 

and the left- correspondence (A_,^_) by 

A_ := -(l-T^i), ^_=T. 

(7 

Note that both are connected by the change a ^ —a. 

The basic polynomial series of the right-correspondence is related with the 
Pochhammer symbols. Its explicit expression is 
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x{x — a){x — 2(7) ... (a; — (n — 1)<t) 



(-a) 




1)! 



1)! 



(to 



if TO < 

m > 
m < n 
TO > 

TO > n 



(11) 



where x — ma with to G Z. This series of polynomials has some interesting 
properties. First of all we can see that the n-th order zero of x" splits into n 
first order zeros of x^""^ at the right of the origin. If the right-correspondence 
is applied to a Taylor power series then, for the positive points of the domain, 
it will be truncated and therefore will converge, but not for the negative ones, 
and therefore could diverge. 

The left-correspondence {A„,^„} is nearly similar to the right one. The 
basic sequence is related with the rising factorial (instead of the falling factorial 
of the previous case), and it is mirror symmetric respect to the y-axis to the 
right-correspondence: in fact, the negative points of a Taylor series discretized 
in this way, always converge. 

According to the formula pop there is only one symmetric representation of 
second order on T (5'-representation) determined by the pair {As,^^}, where 



A. 



1 

2^ 



13s = 2{T + T- 
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and its basic sequence is 



= x[x - [n - 2)a][x - {n - A)a] ...[x+{n- A)a][x + (n - 2)cr] 



_ ri-2 

= x Jl [x + (2i - (n - 2))ct] = a"™ J| [m + 2i - (n - 2)] 



i=0 



1=0 



(sign(m)(T)"|m| 



2)! 



(|m| -n)!! 



(_l)^("-|™|-i)(sign(m)cr)"m 
x(|m| +n- 2)!!(n- |m| - 2)!! 



if n < |m| 
if 
if 



n > |to| 
p(n — m) 
n > |m| 
p(n — m) 



_1 (12) 



where p(n) = (—1)" is the parity function and the double factorial is defined as 



n!! 



n{n — 2)!! with n!! = 1 when n is or 1. 



The behaviour of this basic sequence is much more complicated than those of 
the right and left-correspondences. The zero at the origin of order n of the con- 
tinuous power function of degree n splits into n simple zeros after discretization. 
These zeros appear at the points x = ma with the conditions 



p(n — m) 



-1. 



n > \m\. 



(13) 



Thus, applying the S'-correspondence to a Taylor series we find that at the points 
obeying the conditions (fTS]) there is a cut-off in the series and, hence, the sums 
will always converge in case the function has a well defined parity. 

In the following diagram we display the basis elements (^", A, x*-"^) of the 
three usual representations: 



Right representation Left representation Symmetric representation 

= (XT-i)" ^ ^'l = {XTf ^ ^ (2X {T + T-^y^y 

A+ = i(r-l) ^ A_ = i(l-T-i) ^ As = ^{T-T-^) 

n— 1 n— 1 n—2 

x^+'> ^ l[{x - la) ^ x^^'> = l[{x + la) ^xi"^ = xjjix + {2i - n + 2)a). 

i=0 i=0 i=0 
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19. 







□/ 
/ 



Q- 



J3' 



(a) Some powers for the position in the S-correspondence 



o o 



□ /O,,''' 



/ □ 



□ / 
/ 

/ 

/ 

□ / 

(b) Some powers for the position in the left- (c) Some powers for the position in the 



correspondence 



right-correspondence 



Fi gurG 1 '. Represented in dashed hne the continuous powers and in dots the powers for the 
discrete position with red circles for and blue squares for . The number of zeros is equal 
to the power n for the right and the left cases. 

5 The discretized exponential function and the 
trigonometric functions 

We start with the exponential function which in some sense, is the simplest case. 
The continuous exponential function, as it is well known, is defined through the 
following differential equation 



-r- Vk = kyu, yk{0) = 1, k e 
ax 

and its power expansion series around the origin is 



71=0 

As we mentioned before, the umbral correspondence gives rise to the same 
function either by solving the discrete difference equation, obtained by applying 
the umbral correspondence to the continuous differential equation, or by ap- 
plying the umbral correspondence to its solutions as long as their power series 
converge. Hence, we have to calculate 



(14) 



n=0 



8 



If we use the right-correspondence pT|) the umbral-exponential becomes 

m! 



exp^(fc, ctto) 



ji=0 



n\{m — rt)! 



if m > 



(— m + n — 1)! 



m — 1)! 



if m < 



= (1 + kaY 



For the left-correspondence we get 



exp_(fc, ctto) = < 



ri=0 



n=0 



n!(m — 1)! 

, (-™)! 
n!(— m — n)l 



if m < 



In Fig. 2(a) it is displayed the continuous exponential function together with its 
discrete version by the right-correspondence. It is worthy to note that if in the 
continuous case we found that exp(A: x) — exp{{—k)(—x)), now in the discrete 
case we have that exp_|_(fc,x) = exp_(— fc, — a;). 

Another possibility is to apply the S'-correspondence (fT2|) . So, we obtain 



(ka) 



Yl (m + 2i- {n - 2)) = (ka + {kaY + 1 



4=0 



These discrete exponentials only converge when the continuous variable k 
obeys, [kaY < 1, for all the three correspondences studied in the previous 
sections. 

Also, we can easily study the trigonometric functions (as well as the hyper- 
bolic functions) taking into account their expressions in terms of the exponential 
functions. In this way, we can obtain the discrete trigonometric functions from 
the discrete exponential function. We display in Fig. |2(b)] the sine function. 

We can see that the symmetric discrete trigonometric functions maintain the 
periodic behaviour of the continuous ones. However the period is different 
to the continuous hmction A = 27r/fc. We find it easier in the case of the sine 
function through sins(fc, la) = sins(fc, Oct). We have fixed X ^ la with / G Z, but 
it can be generalized without problems to Z e M. 

sin,(fc,ZCT) = ^ N - [kaY +ika'j' - (^^/l - (kaY - ika'jj = = sm,{k,Oa) 



=> - {kaY + ika'^ = (^^/T^^{kaY - ika 
to solve this identity we define 



^/l - {kaY + '^ka = re'^ 



and this guide us to the solution 



arctan 
1 



ka 



sm{0l) = ^ 6i„ = 



1 mr 
=> kn = - sm — . 

CT / 
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As there are two zeros in each period of the function, we can find the relation be- 
tween the number of waves and the wave length. We have added an superscript 
to remark that these relations are different for each representation, then 

/c'* = _sm — , A'' = — -. (15) 

It is easy to demonstrate the periodic behaviour of the symmetric trigonometric 
functions, in fact introducing (jlSp into de symmetric sinus 

■ r,s ^ . / 27r 
sins(fe , mcr) = sm 1 ™~ 

For the right and left correspondences the trigonometric functions are not 
periodic. However the extremums of the functions and the intersection points 
with the abscissa axis appear periodically and then we can define a wave function 
through this property. With a similar computation than in the symmetric case, 
we find 

1 /27r\ 27rcr , , 

fc+ = - tan — , A+ = — — . 16 

a \ ' / arctan(fc+cr) 

The wave length is bigger in the right and left discretizations and smaller 
in the symmetric discretization than in the continuous trigonometric functions. 
Also in these three correspondences the trigonometric functions have a minimum 
period, it appears in the convergence limit, when ka = I. For the right and the 
left discretizations we have A^^j^ = Scr, in fact it is a 8 point wave, and for the 
symmetric one A^j„ = 4cr, in fact it is a 4 point wave. 

We set X = la and without loss of generality we fix I integer. Then it is easy 
to calculate the change in the amplitude in the right and left correspondences 
through (HH) 

1 ( ^tt\ ( O-n 

sin+(fc+,m(T) = - ((l + ifc+ff)" - (1 -ifc+d)") =cos~" ( 
and 

sin_|_(fc, ma)An(V) = sin+(fc, (rn + nl)a) 
Hence, the non periodic part is then the multiplicative factor 

, /27r 



A„(0 = [see [jjj ■ (17) 

This is also true for the cosine function and for any / G M. Although we have 
a wave length X — la for the trigonometric functions in the right and left 
representations, we can see that they are not periodic. 

It is important to note that all the properties that obey the exponential 
and the trigonometric functions are maintained after the discretization with an 
exception: any property that includes two or more different constants k, because 
they give similar effects than a different lattice spacing a, nor a different x like 
in the continuous case. Thus, for example 

&cpu{k, ma) expij{k, na) = exp[;(fc, (to + n)a), 

exp^(A:, ma) exp[;(fc', ma) ^ expu^k + k' , ma). 



10 



(a) Exponential function with a = 
0.2. 




(b) Sinus function with cr = 0.3 

Figure 2: We represent here the continuous functions with a black solid line, and the discrete 
ones with red triangles the right correspondences and with blue circles the symmetric ones. 



6 Umbral Schrodinger equation 

Let us go to discretize the Schrodinger equation with a potential V{x). We will 
profit the well known fact that if the potential is independent of the time then 
the continuous Schrodinger equation is separable. Thus, from the Schrodinger 
equation 

ihdt(l){x, t) = -—dl(t){x, t) + V{x)(l,{x, t), 
Zm 

and considering (f){x, t) — ip(x)(p{t) we obtain two uncoupled ordinary differential 
equations 

+ V{x)^{x) = Ei:{x), 

ihdtifiit) = Eip{t), 
that we can rewrite with natural units as 

dlipix) + V{x)^{x) = E^Pix), 

(18) 

-idtifiit) = Eifiit). 
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The solution of the equation (|18b ) is 



^{t) = Ce'^\ (19) 

6.1 Constant potential 

In relation with the first equation of psp . the simplest case appears when the 
potential is constant (V{x) = Vb). The solution is 

i^ix) = Ae'*^^ + Be-"''^ if E > Vq, 

(20) 

V'(x) = Ae'^^ + Be-'==" if E <Vq, 



with fc = y^\E — Vol- Therefore, the solution </)(a;, t) is a linear combination 
of trigonometric or exponential functions and, hence, their discretizations have 
been studied in section O 

The energy, for a fixed momentum fc, is the same than in the continuous 
case for any discretization as we can see: 



^ (i/c)" 



H^p = (-A^ + yo)2^^(^ + (-i)"i?)r 



n=0 



n=0 



In spite of this result, the relation between the momentum and the wave length 
X — la is different in the discrete case: 

27r , 1 /27r\ , 1 . /27r 
k = — , fc± = - tan — , /Cs - sm — 
A a \ I J (7 \ I 

Then the phase velocity changes, although the group velocity is the same as in 
the continuous case and, for this reason, the relation between the energy and 
the velocity. 

Let be i = nr where n G Z the discrete time, and r the fundamental time 
length. Following the limits of convergence of the discrete exponentials we find 
two upper limits for the energy, one in relation with the temporal part of the 
wave function (|19p and the other one in relation with the spatial part of the 
wave function ([20]) . 

In fact, we have to take the minor of them in each situation. Taking a » 
Ipianck — 1-62 10~^^m and r « tpianck = 5.39 10~'**s we find values of the 
maximum energy for the time wave Emaxt f-nd for the spatial wave for a proton 
Ep+^max, and an electron E^- ^,nax, 

Emaxt ~ Epianck = 1-22 lO^^eV , E^ax, = 1-38 10^° — cV 

m 

Ep+,„,ax, = 7.94 lO^^ey , E,-^^,,^ = 1.46 lO'^^eV. 
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Remark that these are non-relativistic effects. We can study also the change in 
the ampUtude for the right and left correspondences. The factor of expansion 
or contraction is 



Anil) = 




where I is the number of points into the wave length and n is the number of 
wave lengths along the wave. For the temporal part of the wave function of a 
particle which travel into the LEP, that is during 9 10~5s with an energy of 
1.99 we find that the factor is An{l) - 10^^\ For the spatial part of 

the same particle, an electron, the factor is of order An{l) ^ 3. These factors 
grow a lot with the energy or the mass of the particle. As at this energy the 
relativistic effects are important and we are using now non relativistic quantum 
mechanics, it is comprehensible to find these too big numbers. 



6.2 Infinite potential well 

A little more complicated is the infinite potential well, although it is very similar 
to a constant potential, it reduces the space to a finite region and due to this 
fact new differences will appear between the continuous and the discrete case. 
The potential is 

f if a; e [0, L] 
\ oo if a; ^ [0, L] 



Vix) 



The wave functions must obey the boundary conditions 4^(0) — i^{L) = 
which lead to a countable number of solutions: 

/ a; \ nil 
Wn\x) = sm (^"-^^J 1 ^ ^ ~^ 1 71 e Z, 

in the same way we apply the boundary conditions to all our discretizations. We 
can find the quantum rule for both the right and left correspondences following 
the computation of ([TS]). So 

,1 n 
kZ = — tauTT — 
" cr M 

where L = Ma. Therefore there are only M/2 possible states instead of the 
infinite states of the continuous case, with a degeneration relative to the parity. 
The energy is En = fc^ and it is bounded at n = (M — l)/2. So the maximum 
energy is 



Af - 1 



2M 



At n = M/2 we find a non physical state of infinite energy and infinite wave 
function amplitude. We can do a small computation finding that an electron 
in an infinite well of wide the Bohr radius {AI ^ 10^^) will have a maximum 
state energy of E^^^ ^ lO'^'^eV. The maximum energy grows when increasing 
the number of points M and also when decreasing the mass of the particle. But 
the wave function, as we saw previously, has a convergence range that limits the 
number of states to n e (— M/4, M/4) because k < 1/<t. Hence, the maximum 
energy is much more smaller. 
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We have that the quantum rule for the iS-correspondence is 



,s 1 • " 
fc„ = — sinTT— . 

" cr M 

Again we have a finite number of states, M /2, rounding down, with a degener- 
ation relative to the parity also. The energy is always bounded by h/(2ma'^), 
just as we found with the convergence condition 

In Figure [3] are drawn the energy levels and the wave functions of the ground 
state and the two first excited states for the continuous case and the three 
discretizations. We can see the great fit of the symmetric one for the first 
excited states, but remarking that as its number of states is finite, M/2, then 
upper these wave functions there are not more states. The right and the left 
ones have an asymmetry to the left and right directions of the space and then 
are much more different to the continuous wave functions. 




Figure 3: Energy levels and wave functions of the infinite well. The energy 
levels are in the left, ones in solid line for the continuous case and dashed for 
the discrete case, the shorter red lines are the right/left case energies and the 
longer blue lines are the symmetric case energies. Wave functions are in the 
right, in black solid line are for the continuous case, blue circled points are for 
the S-correspondence, red triangled ones for the right correspondence and green 
squared ones the left correspondence. 
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7 Conclusions 



In this paper we have shown how to apply the techniques of the Unibral Calculus 
for the discretization of equations preserving the point symmetries. We can see 
the advantage of a tool that may be used not only on the equations, but also 
over their solutions. Due to this fact we have not to solve difference equations 
because we alternatively can sum series. Moreover, in the most difficult cases 
we can use numerical methods to get a solution. 

In the umbral discretization we can choose between infinite umbral rep- 
resentations. However, still in the simplest cases they present very different 
behaviours. We can separate all of them in two classes, one preserving the par- 
ity of the equations and solutions (like the S'-representation) and the other one 
that does not (like the right and left representations). From the point of view 
of the local symmetries in Ref ([2]-[n]j [IS]) it was proved that independently of 
the correspondence we get the same discrete point symmetries. 

The study of the discrete Schrodinger equation together with other potentials 
(Coulomb, harmonic oscillator, . . . ) are in progress [31j-|32j and the results will 
be published elsewhere. 
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